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Abstract 

The multi-symplectic form for Hamiltonian PDEs leads to a general 
framework for geometric numerical schemes that preserve a discrete 
version of the conservation of symplecticity. The cases for systems or 
PDEs with dissipation terms has never been extended. In this pa- 
per, we suggest a new extension for generalizing the multi-symplectic 
form for Hamiltonian systems to systems with dissipation which never 
have remarkable energy and momentum conservation properties. The 
central idea is that the PDEs is of a first-order type that has a sym- 
plectic structure depended explicitly on time variable, and decomposed 
into distinct components representing space and time directions. This 
suggest a natural definition of multi-symplectic Birkhoff's equation as 
a multi-symplectic structure from that a multi-symplectic dissipation 
law is constructed. We show that this definition leads to deeper un- 
derstanding relationship between functional principle and PDEs. The 
concept of multi-symplectic integrator is also discussed. 



1 Introduction 

The multi-symplectic Hamilton's Equation is proposed first time by Bridges 

in 1997pp. Birkhoff's Equation is named by Santilli|2I, who just discussed 1- 

dimensional cases. We return to multi-symplectic Hamilton's Equation, and 

we just discuss 1 -|- 1-dimensional cases. It is a local concept of symplecticity 
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for Hamiltonian systems. The concept can be abstracted as follows. Let M 
and K be any skew symmetric matrices on W^[n > 3) and Let S : M" — > 
R be any smooth function. Then we will call a system of the following form 
a Hamiltonian system on a multi-symplectic structure, 

Mzt + Kz^ = y,S{z,x,t), zGE", (x,^)GM^ (1) 

where the general gradient is defined with respect to the standard inner 
product on M", denoted by < •, • >. Associated M and K are the so-called 
pre-symplectic forms 

uj{U,V) =< MU,V >, k(U,V) =< KU,V > . (2) 

The two 2-forms uj and k define two symplectic structure on M™' (m = 
rankM < n) and {k = ranki^ < n) respectively. The pair {uj,k) 
satisfies conservation of symplecticity 

A multi-symplectic integrator is numerical approximation of equation Q 

Mdf'^Zij + Kdx'^'^Zij = {\/zS{zij,Xi,tj))ij, where ztj = z{xi,tj), (4) 

which has a discrete conservation law as 

9f '^Wi J + J = 0. (5) 

The problem with this type to abstract and to discrete the systems 
or PDEs is that it is subject to conservative systems or nonconservative 
systems without dissipation terms. A more general multi-symplectic struc- 
ture is to represent some nonconservative systems with dissipation terms, 
since, — different from form — , symplecticity may depend on time vari- 
ables explicitly, and the right hand is a gradient function (vector function) 
plus a dissipation terms. 

Similarly, let M{z,x^t) and K{z,x,t) are skew symmetric metrics on 
> 2, and B; M" xM^ — be any smooth function. Then a Birkhoffian 
system under the following multi-symplectic structure is given 

dF dG 

M{z, X, t)zt + K{z, X, t)z^ = S7zB{z, ^' *) + ^ + (6) 
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where M and F, K and G are independent. We will show their relationship 
in the next section. 

Structure © is interesting for two reasons: it is a natural generalization 
of structure as well as Birkhoff 's representation for ODEs to PDEs, and it 
can be used to develop symplectic numerical integrators for nonconservative 
systems with dissipation terms. For example, consider the equation describe 
a linear damped string 

utt -Uxx + u + aut + = 0. (7) 

we introduce the new variables p and q, defined below, and equation Q is 
rewritten as first-order form 

ut = p, Ux = q, pt-qx + u + 2p = 0, (8) 

or, with z = {u,p,q)^ € M^, we obtain a multi-symplectic Birkhoff 's equa- 
tions (dsl.6), where 

e°*-^^ 0\ / -e°*-/3^^ 

M = I -e°*-/3^ , = I , (9) 



and 



0/ \e"t-P^ 



i7=(_ie"t-/3-'p,ie-*-/3-^,0r,0), 
G = (^e"*-^^g, 0, ie"*-^^u, 0)^), (10) 
B = -hi^^'l^^[u^ +p2_q2^ + puq]. 

We now show an outline of the paper. In section 2 we introduce ac- 
curate definition of multi-symplectic Birkhoff 's equations. From the point 
of inverse variational problem we discuss the conditions of self-adjointness 
for Birkhoff 's equations. This leads naturally to a variational principle. In 
section 3, we introduce continuous and discrete dissipation laws of sym- 
plecticity, energy and momentum, then give definition of multi-symplectic 
Birkhoffian integrator. In section 4, we consider equation ((T)) to illustrate how 
a multi-symplectic Birkhoffian integrator to preserve the discrete version of 
the dissipation laws. 
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2 Multi-symplectic Birkhoff 's Equation 



The general multi-symplectic structure we shall study is given by 

dF^{z,x,t) _ dFfj_{z,x,t) dzu r dG^{z,x,t) _ dGfj,{z,x,t) dz^, 
^ dzt" dz" ^ dt ^ dzt" dz" ^ dx 

_ dB{z,x,t) dFf,{z,x,t) dG^{z,x,t) ^ 

^ dz^ dt dx ^ ~ ■ 

Equations (jlip are called multi-symplectic Birkhoff 's equations. Why using 
this name is reviewed in Ref.p]. The function B{z, x, t) is called Birkhoffian. 
A representation of a system via multi-symplectic Birkhoff 's equations is 
called a multi-symplectic Birkhoffian representation, and the system is called 
multi-symplectic system. 

Now we consider a general first-order form in the self-explanatory unified 
notation 

M^^{z, X, t)z'^ + K^^z^ + Dmu{z, X, t) = 0. (12) 

We give the conditions of variational self-adjointness for equation (|12() as 
follows. 

Theorem 2.1. (self-adjointness of the covariant general form) A necessary 
and sufficient condition for equations to be self-adjoint in a region 3? 
of points (z,x,t) is that all the following conditions 

iVV + Mvt^ = 0> + K^^ = 0, 

di^ dK^ dK^ dK^ 

dz^ dzi' dz^ ' dz-^ dzi^ dz-" ' (^3) 
dM^^ _ dD dD dKf,^ _ dD dD 

dt ~ dz^ ~ dzi"' dt ~ 'dzy~ Si^*' 

are identically verified in 

Note that the conditions of self-adjointness are formulated for the co- 
variant form because this is the form derivable from a variational principle 
which is the most general possible first-order action functional given by the 
pfaffian action 

£{E) = j j dtdx[Fv{z,x,t)z'i + Gy{z,x,t)z''^- B{z,x,t)]{E). (14) 
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Strip away the integral from £{E) in (|14|) . define a general Lagrange 
density C = Fu{z, x, t)z^ + G^iz, x, t)z^ — B[z, x, t). It doesn't exist a Leg- 
endre transformation (including conjugate momenta ) leads reduction from 
the Lagrange equation to Hamilton's equation. While from the variational 
principle we can write a Birkhoff 's equation directly, and for any self-adjoint 
covariant system, a variational principle can be constructed. 

Proposition 2.1. (Self-adjointness of Birkhoff's Equations). Necessary 
and sufficient condition for a general first order system (ds2.2) to be self- 
adjoint is that it is of the Birkhoffian type, i.e., 

M^i,{z, X, t)z1 + Kfj_yZ'^ + Dmu{z, X, t) = 

dFy{z,x,t) _ dFf^{z,x,t) dzt, dGy{z,x,t) _ dGfj,{z,x,t) dzi, 
\- Q^t, Q^u ^ dt ^ dzt' dz'' ^ dx (^^) 

_ dB{z,x,t) dFf,{z,x,t) dG^{z,x,t) ^ 
^ dz^ dt dx ^ ~ ■ 

3 Multi-symplecticity Dissipation Laws and Multi- 
symplectic Integrators 

The multi-symplectic Birkhoff's equations determines 2 two-forms uj 
and kappa. Since 

^uj{U, V) =< MUt, V> + < MU, Vt> + < MtU, V >, 
Q (16) 
—^{U, V) =< KU^, V> + <KU,V,> + < K,U, V >, 

one obtains 

— w(?7, V) + g^<U, V) =< MUt + KU^, V > - <U,MVt + KV^> + < MtU + K^U, V > 

=< F,tU + G,^U, V>-<F,tV + G,^V, U > 
= 0. 

(17) 

This statement is true for any solution pair {U, V) and we write is as 

dtuj + d^K = 0. (18) 

This is different from the conservation law ((31), — although they look like 
each other — , since that Mj ^ and 0. Therefore we define this law 
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as a multi-symplectic dissipation law for Birkhoffian system, and in the next 
section this definition wih be verified very natural. 

The numerical framework is start from the Birkhoff 's equation (jllj) with 
the associated dissipation law of symplecticity ()18|). The idea is to develop 
a numerical scheme to approximate equations 1)11(1 and it also preserve a 
discrete version of dissipation of symplecticity. A numerical discretization 
of ()18|) can be written schematically as 

(19) 

~l~ {dxG{zi j,Xi,tj^^ij, — z(xj,tj), 

and dl''' and 9^'"' are discretizations of the derivatives dt and dx respectively. 
A discrete version of the dissipation law can be written as follows 

di'^ujij + = 0, (20) 

where 

LVij =< MijUij,Vij >, Kij =< KijUij,Vij > . (21) 

Definition 3.1. The numerical scheme determined by equation (refds3.5) 
is called a multi-symplectic integrator if discrete law h21\) is a dissipation 
law for equations \2U\) . 



4 example 

We consider the following equation 

utt - Uxx + u + 2ui = 0. (22) 
it has a multi-symplectic Birkhoffian representation that has structure 
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M = I -e^* 0), K = \ Q {) 0|, (23) 

and 



-]^e^\u^ -q^ + 2up]. (24) 



It has a dissipation law 



-^(dpAdii) + -^(dg Adu) = -2dpAdu, (25) 
dt dx 
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using the conditions of self-adjointness, the dissipation law can be rewritten 
as 

/rocddt (e^Mp Adu) + -^ (e^Mg A du) = 0. (26) 

dx 

Approximating the equation H22|) in Birkhoffian structure as fohows 

e''*Pi+l/2,j+l — e~^*Pi+l/2,j Qi+l,j+l/2 — Qi,j+l/2 _ 

It + " -^^+1/2^+1/2, 



Ui+l/2,j+l - Ui+i/2,j |'97\ 
= Pi+l/2,j+l/2, ^^'> 



6t 

^j+lj+l/2 — ^j,j+l/2 _ 

- Qi+l/2,j+l/2, 

where 

Pi+l/2,j+l/2 = \e^\Pi+l,j+l +Pi,j+l) + \(i~^\Pi+l,j +Pi,j), 

■Ui+l/2,j + l/2 = -^{Ui+lJ+l + UiJ+l + Ui+lj + Uij), (28) 

1, , 

Qi+l/2,j+l/2 — -^[Ui+l,j+l + UiJ+l + Ui+ij + Uij). 



Discrete the dissipation law (|23j) directly by using the same scheme, we get 

e'^*dpi+i/2,j+i A dui+i/2j+i - e-^^dpi+i/2j A du^+i/aj 

St ^29) 

dQ'j+ij+i/2 A dnj+ij+i/2 - Mi,j+i/2 A dnj_j+i/2 _ ^ 

From the discrete equations (|^^ . we also can construct the dissipation, 
i.e. (PU]) is a symplectic dissipation law for equations H25() . So the scheme 
determined by equations H25() is a multi-symplectic integrator for Birkhoff 's 
equations of system (|^^ . 
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